Vigorous star formation hidden by dust in a galaxy at z — 1.4 
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Near-infrared surveys have revealed a substantial population of enigmatic 

ON ■ 

faint galaxies with extremely red optical-to-near-infrared colours and with 

^: 

a sky surface density comparable to that of faint quasars [|]]. There are two 
scenarios for these extreme colours: (i) these distant galaxies have formed 
virtually all their stars at very high redshifts and, due to the absence of 

•i-H . 

^ ■ recently formed stars, the colours are extremely red [§] and (ii) these distant 



galaxies contain large amounts of dust, severely reddening the rest-frame 
UV— optical spectrum. HR10 (z = 1.44, Ql], |3|]) is considered the archetype of 
the extremely red galaxies. Here we report the detection of the continuum 
emission from HR10 at 850/im and at 1250/um, demonstrating that HR10 is a 
very dusty galaxy undergoing a major episode of star formation. Our result 
provides a clear example of a high-redshift galaxy where the star formation 
rate inferred from the ultraviolet luminosity would be underestimated by 
a factor up to 1000, and shows that great caution should be used to infer 



the global star formation history of the Universe from optical observations 
only. 

The current knowledge of the global star formation activity in the Universe is mostly 
based on the UV continuum luminosity of high- z galaxies Q. However, a crucial phase 
in the evolution of a galaxy might be a vigorous starburst producing copious amounts 
of luminous stars. In this stage, most of the UV stellar radiation would be absorbed by 
the dust grains and re-emitted in the far- infrared (FIR) , making the galaxy dark at UV 
and optical wavelengths. A cosmic background has been first detected at X Q b s ~200/im- 
2mm || and then interpreted by ||. In addition, a cosmic background at 140/im and 
240/im has been recently found by . The presence of such a background suggests the 
existence of a population of distant dusty galaxies. Several models predict the number 
of dusty galaxies at high-z [Q, ||, |H], Q . However, the presently known dusty galaxies 
at z > 1 are limited to a handful of extremely rare objects such as radio galaxies, 
quasars and JiL4S'-selected luminous galaxies ( fill , anc ^ references therein), and 
the existence of a general field population of dusty galaxies has not been established 
yet. Recent surveys for field galaxies found faint (R >24) and very red objects with 
colours typically in the range of 7< R — K <8 [[[], [0|. Their sky surface density, 
~0. 01-0. 02 arcmin -2 at K < 20 [|T], [L3|, is comparable to that of faint quasars [1|], but 
about one order of magnitude less than that of field galaxies with z > 3 fl4 |. 



In order to test if a fraction of these galaxies have these extremely red colours because 
of strong dust reddening, we started a programme aimed at detecting their (sub)-mm 
dust thermal emission. Our first target was HR10, one of the reddest galaxies known 
to date, with / = 24.9, R — K ~ 8, I — K = 6.5 [l[ and the only with a measured 
redshift (z = 1.44 ||). Its UV-optical spectral energy distribution (SED) suggested the 
presence of dust reddening |T^] . HR10 was observed and detected with the Institut 
de Radioastronomie Millimetrique (IRAAQ 30m antenna at 1250/im as well as with 



the James Clerk Maxwell Telescope (JCMT)^ 15m telescope at 850/im in March and 
September 1997 respectively. The details of the observations and data reduction can 
be found in the legend of Figure 1. Throughout the paper we assume Hq = 50 kms -1 
Mpc^" 1 , q = 0.5, and we define h 50 = H o /50. 

Figure 1 shows the SED of HR10. Synchrotron radiation emission as the source of 



the sub-mm fluxes can be ruled out because the radio fluxes |3], [16] constrain the 
extrapolated synchrotron flux at 1250/im to be <0.34 mJy, more than an order of 
magnitude less than the observed one. We therefore interpret the continuum emission 
observed at 850/zm and 1250/im as due to thermal dust emission. Using the ISO upper 
limit at 175/im [I7||, and for dust emissivity indices (3 — 1 — 2 [|Tj|, the temperature 



is found to be in a range 18< Td <45 K. These temperatures fall into the range of 
those of active and star-forming galaxies ( fli~9"l, 12|, [TT| and references therein). The 



total dust mass can be estimated as p0fl : Md = S Vobs D 2 L /[{l + z)Kd{y Tes ^)B{y Test , Td)}, 



where v ohs and v Test are, respectively, the observed and rest-frame frequencies, S Uobs is 
the observed flux density at z/ b s , Dl is the luminosity distance, B is the black-body 
Planck function, T d is the dust temperature and n d = 0.67(^^/250 GHz)' 3 cm 2 g -1 



is the adopted mass absorption coefficient In the two extreme cases considered 

in Figure 1, the total dust mass is M d = 7.3 x 10 8 M and M d = 3.3 x 10 9 
M Q for Td=45 K and T^=18 K respectively. The dust mass increases or decreases of 
a factor of about 2 if (3=1 or (3=2 respectively. If we assume M^ 2 /M^=100 pl |, the 
mass of molecular hydrogen is of the order of M H2 w 10 11 M . The total FIR 
luminosity is estimated by integrating the grey-body curves in the range X rest =10- 
2000/im, LpiR^rest = 47rD|(l + z) 2 J S^obsdu, where -Dl is the luminosity distance and 
S u ,obs is the flux density. We find that L FIR>rest = 3.8 x 10 12 L Q and L FIRjrest = 

1 The JCMT is operated by the Joint Astronomy Centre on behalf of the Particle Physics and 
Astronomy Research Council of the United Kingdom, the Netherlands Organisation for Scientific 
Research, and the National Research Council of Canada 
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1.5 x 10" hw L© for T d =18 K and T d =45 K respectively. 

At low redshifts, HR10 is comparable to the ultra-luminous infrared galaxies (ULIGs) 
which are dusty systems with Lpm > 10 12 L and typically T d ~ 30-60 K [|12|], but with 
dust masses M d ~ 10 7-8 M m lower than those inferred for HR10 ([1121 and references 



l d ~ 1U IVlQ 

therein). At high redshifts, HR10 can be compared to 77^4510214+4724, an ULIG at 
z = 2.286 [2^] whose properties are severely affected by gravitational lensing. If the 
correction for lensing magnification is taken into account (~ 30 x in the infrared p3||). 
77L4S10214+4724 has L FIR = 4.4 xlO 12 L and M d = 1.1 xlO 8 M comparable 
to those of HR10, but it has a warmer dust temperature (40 < T d < 80 K, [|l|, |i~lfl 
and references therein). The warmer dust could be due to the additional source of UV 
radiation provided by the hidden quasar present in the nucleus of 77L4S10214+4724 



24 . 



If we apply the relation between the star formation rate and the FIR luminosity : 
SFR = mQ- 10 L FIR /L & M & yr~\ where * = 0.8 - 2.1 ([f25| and references therein), 
and we adopt \& = 1.5, we derive SFR ~ 570 — 2250 h^Q M yr _1 , which is larger than 
the typical SFRs of \ow-z ULIGs (~10-100 M ?/r _1 []12|), but comparable to that of 



7iL4S10214+4724 (^660 M^r' 1 ]TTJ). If HR10 contains an obscured AGN heating 
the dust, then only a fraction of Lpj R can be ascribed to star formation. The Ha 
line detected in a very low resolution spectrum of HR10 0] has a signal-to-noise ratio 
insufficient to establish the presence of a broad component of this line, which would 
imply that the nucleus of the quasar is not severely obscured. However, the J^-band 
morphology of HR10 is spiral-like and not strongly nucleated, and it does not support 
the presence of a directly visible quasar nucleus. 

Due to the steep rise of the grey-body dust spectra towards X rest ~ 100 — 200/im, 
the observations at longer wavelengths make the sub-mm and mm flux (at a fixed 



observed A) of a dusty galaxy roughly constant for approximately 1 < z < 10 ||, 
|T[| . For instance, we moved HR10 to three different redshifts (z = 0.7,2.2,3.3) and 
we computed how its properties would change accordingly adopting /3 — 1.5 and an 
optically thin grey-body : S u oc b ,l3 B u (T). We find that the inferred dust temperature 
would increase from ?al5-20 K at z = 0.7 to ~ 45-50 K at z = 3.3. Lpm and Md would 
be ~ 0.1, 7, 25 x 10 13 L Q and ~ 30, 10, 5.6 x 10 8 M for z = 0.7 (T d = 18 K), 
z = 2.2 (Td = 35 K), and z = 3.3 (Td = 45 K) respectively. The main implication is 
that a galaxy with observed (sub)-mm properties as those of HR10 would be classified 
a dusty ultra-luminous infrared galaxy in a wide range of redshifts. 

The discovery of dust associated with HR10 has several implications. First of all, it 
sheds light on the nature of this galaxy, demonstrating that HR10 is a very dusty star- 
forming galaxy where, similarly to low-z ULIGs, most of the energy is emitted in the 
far-infrared. The second implication is generally related to the problem of estimating 
the star formation rate in distant galaxies. The global history of the star formation in 
the Universe is presently inferred mostly by the continuum UV luminosity of optically- 
selected high-z galaxies [|J, from which the derived SFRs are in the range of 4-25 
h^Q M Q yr _1 [TJj]. However, such estimates can be severely hampered by the presence 



of dust extinction which modifies the shape and reduces the flux of the UV spectra. 
Observations of low-redshift galaxies have also pointed out the limits of the UV-blue 
light as a SFR estimator (see for instance [EB| and references therein). HR10 allow us 



to see an extreme example of this problem applied to a high-z system. In fact, the Ha 
luminosity of HR10 would imply a SFR ~ 80h^ M Q yr _1 ||, whereas the luminosity 
of the UV continuum at 2800 A would suggest SFR ~ l/i^~ 2 M Q yr _1 both at least 
one order of magnitude less than the SFR suggested by the FIR luminosity. In this 
regard, our results show that the global star formation history of the Universe can be 

fully traced only if the effects of dust are taken into account. Objects like HR10 would 
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have not been found by neither optical imaging surveys based on the Lyman-continuum 
break [TJJ] or on strong emission lines, nor by IRAS surveys, nor by traditional quasar 



surveys. Instead, our results demonstrate that the combination of optical and near-IR 
deep imaging, coupled with (sub)-mm observations, is an efficient method to find dusty 
galaxies at high-z. Recent deep SCUBA imaging suggests indeed the presence of a 



population of faint sub-mm sources [28 



The observations of HR10 suggest that the star formation in distant objects occurs 
with different modes, and that the most vigorous episodes of star formation probably 



arise in dusty environments as predicted by several models |9], |10[ [27| . However, our 
data cannot tell us if HR10 is forming its first generation of stars, or if the starburst is 
occurring in an already formed system. Nevertheless, in both cases we are witnessing a 
major episode of star formation: if the burst lasts about 10 7 years (a typical time scale 
for a starburst Jl2|), the total mass of gas converted into stars is of the order of 10 10 



M , which is approximately 10% of the total mass of a present-day massive galaxy. 

Finally, we find relevant to note that the observed properties of HR10 fit into the 
predictions of || that the sources at z ~ 0.5 — 2.5 which contribute to the cosmic 
FIR — mm background should have fluxes around 10-100 mJy at X b s = 200/im. For 
18 < Td < 45 K, the expected flux of HR10 at X b. s = 200/zm is in the range of 
10-40 mJy. Moreover, we also note that the sky surface density of dusty galaxies at 
X bs = 175/iin predicted by [0] 0.05 arcmin -2 ) is within a factor of five similar to 



that of the extremely red galaxies [fL], [13|]. Future observations will provide clues on 



what fraction of the extremely red galaxies contributes to the FIR — mm background. 
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Figure Caption 



Figure 1: Filled circles: fluxes of HR10 taken from |], ||, [T^, |16| and this work. 



The IRAS 3a upper limits at X b s = 12, 25, 60, 100/im were derived from co-added 

survey data with the SCANPI procedure available at IPAC (Infrared Processing & 

Analysis Center). The curves show two extreme cases of grey-bodies consistent with 

the data. Dashed curve : optically thin grey-body, S u oc v^B v {T^), with (3 = 1.5 and 

Trf=18 K. Continuous curve : grey-body optically thin at v < z/ =750 GHz, S u oc 

B v (T){l-exp[-(u/u G f}} 0, with (3 = 1.5 and T d = 45 K. HR10 was observed for 64 

minutes (on source) with the MPIfR 19-channel bolometer [|9| (1250/xm, 11" (FWHM) 

beam, 32" chop throw, zenith opacity 0.09-0.3). Flux calibration was achieved using 

Uranus as primary calibrator and Mars and pointing quasars as secondary ones. The 

mean atmospheric level was found averaging the measurements of 16 channels and 

subtracted from the central pixel, which targeted HR10. The final flux at 1250/xm is 

4.9 ±0.7 mJy The simultaneous observations at 450 and 850/iin were made with the 

SCUBA bolometer array (Robson et al., in preparation) at the JCMT in photometry 

mode (108 minutes on source, zenith opacity 0.3 at 850/im). The sky background 

was removed from the flatfielded and despiked data by subtracting the average signal 

from 5 neighbour bolometers (containing no source signal). The data were calibrated 

using photometry and beam maps of Mars and Uranus. The flux at 850/nn and 450/im 

resulted to be 8.7 ± 1.6 mJy and <180 mJy (3<r) respectively. The la error bars of 

the 850/zm and 1250/xm flux densities indicate the statistical uncertainties, whereas the 

uncertainty on the absolute flux calibration is estimated to be ~15-20% both in the 
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I RAM and SCUBA data. 
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